Crystal structure and properties of barium thorate BaThOs from first principles 
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The phonon spectrum of cubic BaThOs with the perovskite structure is calculated from first 
principles within the density functional theory. The analysis of unstable modes in the phonon 
spectrum enables to determine the symmetry of all possible distorted phases, calculate their energies, 
and show that the ground-state structure of barium thorate is Pbnm. For this structure, the static 
and optical dielectric constants, elastic moduli, heat capacity, Raman spectra, and the energy band 
gap in the LDA and GW approximations are calculated. The possibility of the structural phase 
transitions in BaThOa is also discussed. 
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I. INTRODUCTION 

The ABO3 perovskite oxides constitute an important 
family of materials that display interesting physical prop- 
erties (ferroelectricity, ferroelasticity, superconductivity, 
magnetism) and are widely used in technological applica- 
tions. Barium thorate BaThOa is one of the members of 
this family. High thermal stability of barium thorate (its 
melting point is above 2300°C) and relatively low work 
function of thorium enable to use BaThOs for thermionic 
cathodes of high-intensity discharge lamps (see, for ex- 
ample, patents^"— ). Barium thorate ceramics doped with 
neodymiumii^ and yttrium^ are solid electrolytes, which 
have the highest proton conductivity at temperatures 
above 900°C among solid electrolytes. Since barium tho- 
rate is produced as the fission product in nuclear reactors 
(especially in uranium-thorium reactors) , the knowledge 
of its properties is also important for predicting the reli- 
ability of nuclear fuel elements. 

Insufficient knowledge of the properties of BaThOa be- 
comes evident from the fact that even its crystal struc- 
ture is not well established. In the very first work^ 
the structure of barium thorate was identified as a sim- 
ple perovskite structure with the lattice parameter of 
a « 4.48 A. A more detailed study^i found the superstruc- 
ture reflections on the diffraction patterns, but because of 
the lack of splitting of the strong peaks, the structure was 
considered as a pseudocubic one with a doubled lattice 
parameter. Finally, after the splitting of the main re- 
flections was observed^^ the structure of BaThOa was 
identified as orthorhombic, but neither the space group 
nor the atomic coordinates were determined. 

Thermodynamic properties of barium thorate (the 
Gibbs energy of formation) were determined in Refs. [ll| 
and [l^. In Ref. [l^, the electronic structure, optical and 
elastic properties of BaThOs were calculated from first 
principles using the FP-LAPW approach, but for some 
reason the calculations were limited to the cubic five- 
atom unit cell, whereas X-ray studies clearly indicate the 
lower symmetry. According to these calculations, the cu- 
bic barium thorate is a direct-gap insulator with the band 
gap of Eg = 5.7 eV. This result, however, contradicts 
the density-of-states calculations presented in Ref. [ij. 



according to which the energy gap between the valence 
and conduction bands is 3.33 eV, and the calculations of 
the 62(0-1) optical spectra in which the absorption starts 
at an energy below 5 eV. 

The knowledge of the true crystal structure of BaThOs 
is crucial for correct prediction of its properties. In this 
work, we use the first-principles calculations to deter- 
mine the equilibrium structure of barium thorate, calcu- 
late some of its properties, and discuss the possibility of 
the structural phase transitions in it. 



II. CALCULATION DETAILS 

To predict the properties of BaThOa, we must first de- 
termine its ground-state structure which has the lowest 
energy at T = 0. For this purpose, we must first cal- 
culate the phonon spectrum of its parent Pm3m phase 
and then, by adding the distortions corresponding to un- 
stable modes in the phonon spectrum to the structure, 
seek for a minimum-energy structure in which the energy 
of all optical phonons at all points of the Brillouin zone 
are positive and the structure is mechanically stable (the 
stability criterion is the positive values of the determi- 
nant and all leading principal minors of the 6x6 matrix 
of elastic moduli in the Voigt notation). 

In this work, the first-principles calculations were per- 
formed within the density functional theory (DFT) using 
the ABINIT softwareJ^ The exchange-correlation inter- 
action was described in the local density approximation 
(LDA). Pseudopotentials for Ba and O atoms used in the 
calculations were taken from Ref. [l3- Scalar-relativistic 
pseudopotential for the Th atom was constructed using 
the RRKJ scheme^^ with the OPIUM program.— To im- 
prove the transferability of the pseudopotential, the s lo- 
cal potential correction^^ was used. The parameters used 
for the construction of pseudopotentials are presented in 
Table|Tl Testing of the Th pseudopotential using Th02 as 
an example revealed its high enough quality: the calcu- 
lated lattice parameter of this compound (5.606 A) dif- 
fered from the experimental value by only 0.15% and the 
bulk modulus (205 GPa) differed by 4%. 

The lattice parameters and the equilibrium atomic co- 
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TABLE I. Electronic configuration of atoms and parameters used for construction of pseudopotentials: r^, r^, and are radii 
of the pseudopotential cores for s-, p-, and d-projections; q^, Qp, and qd are the cut-off wave vectors used in the optimization 
procedure; rmin, r-max, and Vioc are the range limits and depth of the correcting local potential. All the values are in Hartree 
atomic units except for the correction energy which is in Ry. 



Atom 


Configuration 








Is 


qp 


qd 


^min 


^max 




Ba 




1.85 


1.78 


1.83 


7.07 


7.07 


7.07 


0.1 


1.68 


1.95 


Th 




1.88 


2.04 


2.04 


7.57 


7.27 


7.07 


0.01 


1.76 


0.75 


O 


2s^2p^?,d° 


1.40 


1.55 


1.40 


7.07 


7.57 


7.07 









ordinates in the unit cell were determined from the con- 
dition that the residual forces acting on the atoms are 
less than 5 • 10"^ Ha/Bohr (0.25 meV/A) and the total 
energy is calculated self-consistently to an accuracy of 
better than 10~^° Ha. The integration over the Brillouin 
zone was performed on the 8x8x8 Monkhorst-Pack mesh 
for the cubic phase or on the meshes with equivalent fc- 
point density for low-symmetry phases. The maximum 
plane-wave energy was 30 Ha. 

The quasiparticle band gap of BaThOa was calculated 
using the so-called one-shot GW approximation! 
The Kohn-Sham wave functions and energies calculated 
within DFT-LDA were used as a zeroth-order approx- 
imation. The dielectric matrix ecc (q, was com- 
puted for a 6x6x6 q-point mesh from the independent- 
particle polarizability matrix Pqq, (q,a;) calculated for 
3743 reciprocal-lattice vectors G(G'), 20 occupied and 
280 unoccupied bands. The dynamic screening was 
described using the Godby-Needs plasmon-pole model. 
The components of wave functions with kinetic energy 
up to 24 Ha were used in these calculations. The en- 
ergy correction to the DFT-LDA solution was computed 
as diagonal matrix elements of the [E — E^d operator, 
where S = GW is the self-energy operator, E^c is the 
exchange-correlation energy operator, G is the Green's 
function, and W — e~^v is the screened Coulomb in- 
teraction. In the calculations of E, the components of 
wave functions with kinetic energy up to 24 Ha for both 
exchange and correlation parts of S were used. The ac- 
curacy of the band gap calculation estimated from the 
convergence tests is ~0.05 eV. 



III. RESULTS 

The phonon spectrum of BaThOa in the cubic Pmim 
phase is shown in Fig. [TJ It is seen that there are nine 
unstable modes in the phonon spectrum of barium tho- 
rate. They include the triply degenerate Fis mode as- 
sociated with the ferroelectric instability; the triply de- 
generate mode describing the uniform deformation 
of the oxygen octahedra; the triply degenerate i?25 mode 
and nondegenerate M3 mode, which describe the rotation 
of the octahedra; the doubly degenerate antiferroelectric 
X5, Xg, and M5 modes; the non-degenerate mode 
describing spatially non-uniform deformation of the oc- 
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FIG. 1. Phonon spectrum of BaThOs in the cubic Pm3m 
phase. The labels indicate the symmetry of unstable modes. 



tahedra; and the triply degenerate Ri^ mode describing 
the vibrations of the Ba atom in its oxygen environment. 

The energies of all phases obtained from the conden- 
sation of the unstable phonons as well as the energies 
of Pbnm and Cmcm phases resulting from simultane- 
ous condensation of the R25 and M3 modes are given in 
Table HIl Among these phases, the Pbnm phase has the 
lowest energy. 

The calculation of the phonon spectrum in the Pbnm 
phase shows that the frequencies of all optical phonons 
at the center of the Brillouin zone and at the high- 
symmetry points on its boundary are positive; the de- 
terminant and all leading principal minors constructed 
from the clastic tensor are also positive. This means 
that the Pbnm phase is the ground-state structure of 
BaThOs. The calculated lattice parameters of this phase 
are a = 6.3140, b = 6.4039, and c = 8.9578 A; the atomic 
coordinates are given in Table [HTI The obtained lattice 
parameters are in good agreement with the experimental 
data of Ref. i (a = 6.345 ± 0.002^^5 = 6.376 ± 0.002, 
c = 8.992 ± 0.002 A) and Ref. M {a ^ 6.35 ± 0.01, 
6 = 6.387 ± 0.009, c = 8.995 ± 0.009 A). 

We now consider some physical properties of BaThOs 
in the ground state. Although the ferroelectric instabil- 
ity was found in the parent cubic phase, the structural 
distortions suppress this instability, and in the ground 
state barium thorate is a paraelectric. The static di- 
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TABLE II. The relative energies of the low-symmetry phases 
of BaThOs resulting from the condensation of unstable 
modes. The most stable phase is in boldface. 



Phase 


Unstable mode 


Energy (me 


Pm3m 


— 





PA2/mmc 


X3 


-55.9 


Cmcm 




-90.9 


Pmma 


X5 


-99.3 


Pmma 




-101.4 


Cmcm 


Xs 


-118.8 


M/mmm 


R15 


-123.0 


R3m 


Rl5 


-148.0 


P4m2 


T25 


-153.1 


Amm2 




-228.0 


Pmma 




-234.5 


R32 


Fas 


-255.9 


Cmmm 




-303.2 


R3m 




-326.4 


PAmm 


Tis 


-337.2 


PA/mhm 


M3 


-529.2 


I A /mem 


R25 


-599.3 


Cmcm 


R25 + Ms 


-696.4 


R3c 


R25 


-710.1 


Imma 


R25', R15 


-754.4 


Pbnm 


R25 + M3 


-797.8 



TABLE IV. Calculated frequencies of the Raman-active 
modes for the Pbnm phase of barium thorate. 



The mode symmetry 



Frequency (cm ^) 



A, 

Big 
B2g 
B'ig 



75; 89; 157; 207; 300; 341; 523 
95; 109; 169; 283; 312; 387; 666 
92; 217; 331; 508; 681 
96; 123; 298; 520; 613 
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FIG. 2. Calculated Raman spectra for a poly crystalline sam- 
ple of BaThOa at 300 K. 



TABLE III. Calculated atomic coordinates in the orthorhom- 
bic Pbnm phase of barium thorate. 



Atom 



Ba 
Th 
01 

02 



Position 



4c 
46 
4c 
Ad 



y 



-0.01246 
0.50000 
0.10447 
0.70713 



0.03717 
0.00000 
0.46199 
0.29227 



0.25000 
0.00000 
0.25000 
0.05723 



The energies of different phases obtained in this work 
(Table |n| enable to make some remarks on the possible 
phase transitions in barium thorate. The Imma phase 
closest in energy to the Pbnm ground-state structure has 
the energy only by 43.4 meV higher than the ground-state 
energy. This means that a second-order phase transition 
Pbnm Imma can occur in BaThOa at a temperature 
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electric tensor in the Pbnm phase is characterized by 
three diagonal components of — 24.7, e°y = 25.0, and 
= 37.6; the tensor of the optical dielectric constant is 
described by the components of el^ = 4.42, — 4.35, 
and = 4.28. The elastic moduli are Cn = 195 GPa, 
C22 = 200 GPa, C33 = 199 GPa, C12 = 99 GPa, C13 = 
83 GPa, C23 = 71 GPa, C44 = 55 GPa, C55 = 51 GPa, 
Cge = 64 GPa. The bulk modulus in the orthorhombic 
phase is B = 121.9 GPa; in the Pmim phase its value 
is 124.2 GPa and is very close to the value of 124 GPa 
obtained in Ref. |13|. 

For identification of the Pbnm phase of BaTh03 it 
is convenient to use the Raman spectroscopy. The cal- 
culated frequencies of the Raman-active modes in this 
phase are given in Table IIV| and the theoretical Raman 
spectrum calculated for a polycrystalline sample at 300 K 
using the formulas from Ref. |21j is shown in Fig. [51 
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FIG. 3. Comparison of the calculated heat capacity Cv in the 
Pbnm phase of barium thorate (solid line) with the experi- 
mental data22. (points). 
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slightly above 300 K. In the experimental study of the 
heat capacity)^ no anomalies were detected in the tem- 
perature raiige 320-820 K, but the peculiarity of the re- 
sults of Ref. '22 is that the measured values of the specific 
heat are much higher than our calculated ones (Fig. [3]) 
and, more significantly, they even exceed the theoretical 
limit of 3i? per atom in the unit cell (the Dulong-Petit 
law). This may indicate that the investigated sample 
was likely not pure^^ We believe that the predicted phase 
transition was hidden by this fact. In this regard, we note 
that the features in the dielectric properties of BaThOs 
at about 260° C, which were observed in Ref. ^ and in- 
terpreted as an evidence for the ferroelectric phase transi- 
tion, can actually be associated with the Pbnm — > Imma 
phase transition under discussion. 

To resolve the contradictions mentioned in the Intro- 
duction on the energy band gap Eg in cubic BaThOs, we 
performed our own calculations of Eg in both the LDA 
approximation used in this work and in the GW approx- 
imatio n^^i^° which takes into account the many-body ef- 
fects and enables to obtain the energy band gap in much 
better agreement with the experiment According to 
our data, in the LDA approximation E^^^ = 2.907 eV 
in the PmZm phase; in the orthorhombic Pbnm phase 



the band gap increases to 3.947 eV. As the DFT calcu- 
lations always underestimate the energy band gap, more 
realistic Eg values can be obtained in the GW approx- 
imation. These calculations gave the energy band gap 
E^^ w 4.40 eV for the cubic phase. This value is in- 
termediate between the two values, 5.7 eV and 3.33 eV, 
reported for this phase in Ref. In both cubic and or- 
thorhombic phases, the extrema of the valence and con- 
duction bands are located at the F point of the Brillouin 
zone (the optical transitions are direct). 

The calculations presented in this work were performed 
on the laboratory computer cluster (16 cores). 

IV. CONCLUSIONS 

First-principles calculations of the phonon spectra and 
structure of cubic BaThOa and its distorted phases have 
shown that the ground-state structure of barium thorate 
is Pbnm. The static and optical dielectric constants, elas- 
tic moduli, heat capacity, Raman spectra, and the energy 
band gap in the LDA and GW approximations have been 
calculated for this structure. It has been shown that the 
structural phase transition Pbnm — >■ Imma can occur in 
BaThOs at a temperature slightly above 300 K. 
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